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ABSTRACT 

We study a generalization of the Randall-Sundrum mechanism for generating the weak/Planck 
hierarchy, which uses two rather than one warped extra dimension, and which requires no nega¬ 
tive tension branes. A 4-brane with one exponentially large compact dimension plays the role of 
the Planck brane. We investigate the dynamical stability with respect to graviton, graviphoton 
and radion modes. The radion is shown to have a tachyonic instability for certain models of the 
4-brane stress-energy, while it is stable in others, and massless in a special case. If stable, its 
mass is in the milli-eV range, for parameters of the model which solve the hierarchy problem. 
The radion is shown to couple to matter with gravitational strength, so that it is potentially 
detectable by submillimeter-range gravity experiments. The radion mass can be increased us¬ 
ing a bulk scalar field in the manner of Goldberger and Wise, but only to order MeV, due to 
the effect of the large extra dimension. The model predicts a natural scale of 10 13 GeV on the 
4-brane, making it a natural setting for inflation from the ultraviolet brane. 
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1 Introduction 


Brane-world scenarios have undergone considerable theoretical scrutiny in recent times, largely 
due to the novel solutions which they provide for the gauge-hierarchy problem. Three kinds of 
such proposals have been made: the ADD scenario [|IJ, in which two extra dimensions are very 
much larger than ordinary microphysical scales; the intermediate-scale scenario ||, in which 
the scale of gravity and the compactification scale are very close to one another; and the RS 
scenario 0, in which the geometry of the extra dimensions is exponentially warped. 

The RS scenario differs qualitatively from the other two, with its departure from product- 
space geometries potentially implying many new kinds of low-energy features, including possible 
violations of Lorentz-invariance Q, multiple or quazilocalized gravitons [.J, |(J and phenomeno¬ 
logically acceptable modifications of the cosmological Friedmann equations [0, among other 
possibilities. 

Most of the effort on exploring this scenario has focused on 5 spacetime dimensions, making 
it difficult to ascertain which features are generic to warped geometries and which are artifacts 
of five dimensions. Indeed, some features of these models are very likely to be specific to five 
dimensions. For instance, the requirement of negative-tension branes - which are generic in 
five dimensions - does not arise in six dimensions [^]. j9). It is also unlikely that the Friedmann 
equations are modified in six dimensions in the same way as in five. Moreover, the absence of 
Kaluza-Klein excitations of the metric’s radion mode is also a 5-dimensional artifact, due to 
the trivial geometrical nature of 1-D manifolds. 


There have been numerous proposals for higher-dimensional generalizations of the RS idea. 
One of the earliest was to consider intersections of codimension-1 branes as the 3-branes [llO 


Others involved modeling the 3-brane where we are supposed to live, or in some cases where 
gravity is localized, as a cosmic string or higher-dimensional defect 
on warped higher dimensional spaces includes [T2| . 


TlJ. Other relevant work 


A particularly attractive six-dimensional warped model has been considered in various con¬ 
texts by several authors fT3|-[|I(|. This model is related to the AdS soliton |17|, a double 
analytic continuation of a planar AdS Black hole metric, and involves two compact dimensions 
having the topology of a disc with a conical singularity at its centre. The boundary of the 
disc occurs at a (Planck) 4-brane and a (TeV) 3-brane is placed at the conical defect. The 
stress energy of the 4-brane requires an anisotropic form which could arise from the smearing 
of 3-branes around the 4-brane, as suggested in ref. || , or from Casimir energy of light particles 
confined to the 4-brane 
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Since all observable consequences of this (or any other) geometry only involve the theory’s 
low-energy degrees of freedom, essential to understanding its physical implications is a deter¬ 
mination of its low energy spectrum. While this has been partially done in previous references, 
it is the purpose of this paper to provide a complete accounting of the metric modes, especially 
as regards the elusive radion mode. 


We find results which differ interestingly from what obtains for five-dimensional RS models. 
Instead of a massless radion, we find that generically the mass squared is nonzero, and possibly 
negative, depending on details of the 4-brane stress tensor. For a special case involving Casimir 
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energy on the 4-brane, the radion is an exactly massless modulus at the classical level. Quantum 
corrections (which we do not here calculate) might stabilize the radion in this case. If the radion 
mode is stable, the magnitude of the mass squared is of order (10~ 3 eV) 2 . Both the stable and 
loop-stabilized cases might therefore make this mode of interest for table-top tests of gravity. 
The tachyonic instability, if it occurs, does so regardless of the value of the radial size of the 
extra dimensions in the static solution. We show that it is straightforward to cure this problem 
in the manner of Goldberger and Wise |I8|, by adding a bulk scalar field which couples to the 
branes. The mass of the radion, once stabilized, is suppressed relative to the Planck scale by 
an additional fractional power of the warp factor, which puts it in the MeV rather than TeV 
range. Although this could potentially have been problematic, we find that the coupling of the 
radion is similar to that of gravity, because its wave function is not strongly peaked on the TeV 
brane. Therefore, although there are cosmological contraints on this model, it is not ruled out 
by constraints from supernova cooling or radion production in colliders. 


We organize our presentation as follows. In section 2 we will introduce the model at the 
static level. In section 3 we will find the dynamical perturbations for 4-D modes which transform 
as tensors (gravitons), vectors, and scalars (the radion). We will show that, whereas the tensor 
and vector modes have a massless ground state, the radion mass squared is generically nonzero 
and possibly tachyonic, although its magnitude is exponentially small. In section 4 we will 
show how a bulk scalar field can stabilize the radion mode, and discuss the phenomenology of 
the model. A summary is given in section 5. 


2 The AdS Soliton in Randall-Sundrum Models 


In this section we present the AdS soliton |T7| and its key properties relevant for braneworld 
applications. A more detailed description of this spacetime including its role in terms of the 
AdS/CFT correspondence can be found in ref. 0 - 


The (p + 2)-dimensional AdS Soliton first arose in connection with the AdS/CFT corespon- 
dence |L9|] as a double-analytic continuation of a (p + 2)— dimensional planar AdS black hole 
metric. In this context the AdS soliton geometry is relevant for the strong coupling descrip¬ 
tion of Lorentzian signature superconformal field theories in which both supersymmetry and 
conformal invariance have been broken 0. Specifically for p = 3 and p — 5 it provides a 


description of strongly coupled QCD 3 and QCD 4 respectively. 


2.1 The AdS Soliton 

In this article we will be interested in the six dimensional (he., p — 4) AdS soliton for which 
the line element may be written, 

ds 2 = a(r) ( f(r)dd 2 + p^dx^dx + a~ l (r)f~ l (r)dr 2 ( 1 ) 

where the metric functions are given by, 

fir) = j2 ~ and °( r ) = ( 2 ) 
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and r/ /it , is the four dimensional Minkowski metricThis is a solution to six dimensional Einstein 
gravity with negative cosmological constant 


A “ I?~ 5*’ 


( 3 ) 


The spacetime is asymptotically locally AdS as r —> oo; below we will cut off the radial extent 
by inserting a 4-brane at a finite value of r. For convenience we have normalized a(p) = 1, since 
we will be interested in placing the standard model on a 3-brane situated at that position. 


The range of the r coordinate is p < r < oo, with the geometry smoothly ending at r = p 
provided that the 6 coordinate is periodic, with period 


P = 


47 tL 2 


( 4 ) 


This also requires that, in the context of supergravity, fermions are antiperioclic in the 6 di¬ 
rection. It is important to note that this geometry is everywhere smooth and nonsingular 
including at r = p where the circle parameterized by 6 smoothly shrinks to a point and the 
geometry ends. An attractive feature of this geometry is that it ends in a natural and nonsin¬ 
gular fashion, allowing constructions similar to those proposed in refs, 
of uncontrollable and likely unphysical curvature singularities. 


21 , but which are free 


To construct a brane-world model, we will want to imagine that we are living on a (TeV) 
3-brane at r = p, thereby introducing a conical defect there of size 6 = k 2 t 3 , where r 3 is the 
3-brane tension and k 2 is related to the six-dimensional Newton’s constant by Kg = 87 tG 6 . This 
modifies eq. ([4]) to become 

or 2 

(3=-(2k-8). (5) 

5 p 

The extra dimensions are compactified by terminating the space at a 4-brane at r = R. 


In the horospheric coordinate system, the proper distance from p to r along the r-direction 
is given by 


r = 


dr 

VaJ 

and so r/p = cosh(2fcf/5) ~ e r ^ L if r p. 


= k 1 cosh 


-l 


(r/p) 


5 / 2 ' 


( 6 ) 


We will often find it enlightening to express the solution in polar coordinates, r, where the 
line clement has the form 


ds 2 = a(r)j] llu dx fJ, dx 1 ' + b{r)d6 2 + dr 2 (7) 

Here the metric coefficients are given by 

a(f) = cosh 4/5 (AT); 6(f) = b 0 smh ( fcr ) ( 8 ) 

V ' V ' W cosh 6/5 (fcf) V ; 

where bo = k~ 2 if the point f = 0 is regular, and 6 G [0, 27t] . In general we will suppose the 
3-brane has nonvanishing tension located at this point. Then the conical singularity at f = 0 

1 Our metric’s signature is ‘mostly plus’ and we adopt MTW |2^] curvature conventions. 
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introduces the deficit angle given by 2n(ky/b^ — 1). In these coordinates the proper distance 
between two radii f\ and r 2 is simply their difference, f 2 — f \. We denote the radial position of 
the 4-brane by r — R. 


We close this section with some comments regarding the stability of the AdS soliton. Since 
the AdS soliton is constructed from multiple analytic continuations of a black hole space time 
one might worry about dynamical stability of the solution. In general such analytically con¬ 
tinued space times are not always well behaved. For example beginning with the Reissner- 
Nordstrom black hole in asymptotically flat space one can analytically continue the metric 
in such a way as to allow for arbitrarily large negative values for the mass parameter. Solu¬ 
tions such as this are inherently unstable i.e. small perturbations around the background are 
tachyonic—see ref. 0 for a detailed discussion. 


One of the key results of ref. [|lj was that (for p — 3) the AdS soliton was found to be 


perturbatively stable to such linearized fluctuations. Further, in ref. [23| this proof was extended 
to arbitrary p > 2 —see also ref. |24[] for a recent discussion. One consequence of this proof is 


that at least locally within the space of solutions to the Einstein equations with asymptotically 
locally AdS boundary conditions the AdS soliton represents the minimum energy solution. 


It is one of the purposes of this paper to investigate whether the perturbative stability of this 
space time persists when the geometry is cut off by the introduction of the 4-brane discussed 
above. 


2.2 The Gauge Hierarchy 

To understand how this model solves the gauge hierarchy problem, let us imagine that all 
the fundamental scales Mg, k , and 1/R are of order TeV. Then the standard reduction of the 
gravitational action from 6 to 4 dimensions (using polar coordinates) gives the 4-D Planck mass 
as 

M 2 = 

p 


Thus by taking e 3kR,//5 ~ 10 16 , corresponding to kR = 60, we can explain the largeness of the 
Planck scale. 

Notice that the relation M 2 ~ a 3//2 (i?)(TeV) 2 differs from the analogous relation in the 5-D 
RSI model (i.e., the Randall-Sundrum model which is compactified by the presence of the TeV 
brane), M 2 ~ a(i?)(TeV) 2 . The additional factor of a 1 / 2 is coming from the b 1//2 part of the 
measure, which gives the size of the extra compact dimension that was not present in the 5-D 
model. This shows that the present model is a hybrid of the RS and ADD scenarios, in that 
the hierarchy is due to a combination of warping and having a large extra dimension. 

This difference can also be seen by considering the physical mass of a 4-D scalar field which is 
confined to a 3-brane at a fixed position (r, 6) in the bulk. Since we are taking the fundamental 
scale to be TeV, we should assume that its bare mass parameter m is of this order. But the 


Mg / dr dO a(r)Jb(r) 


Mi 


A a 3 / 2 f'R'l — _ 6_ p 6kR/5 


k 2 


\R) = 


Mi 


k 2 


(9) 
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physical mass is determined by the usual argument of canonically normalizing its kinetic term: 


Sk = — - ax a 


a + m 2 (j) 2 


drx 


+ a(r)m 2 (j) 2 


( 10 ) 


Thus the physical mass is given by m 3 = mya(r). If we take the particle all the way to the 

4-brane, its mass does not reach the Planck scale, but rather a smaller one, a(R)~ l ^M p ~ 
10 - 32 / 6 M p ~ 10 13 GeV. This reflects the fact that the strength of gravity is still diluted for a 
4-D observer on the 4-brane, by the large extra dimension. We should refer to it as the “10 13 
GeV brane” rather than the Planck brane. 


2.3 Properties of the Branes 


The 4-brane we need for compactifying the AdS soliton can be constructed by the standard 
cutting and pasting procedure. Here, the metric in eq. (p) will be cut along the surface r = R 
and then pasted onto a mirror image of itself. The resulting space time is then a solution of 


hi 


-2 


V—G ^ Rmn ~ -RGmn — A Gmn^ 


V-gS ab 5 a M s b N 5 


r — R\ 

tfJ’ 


(ii) 


where Gmn is the metric given in eqn.(j^), A = — 10 L 2 is the cosmological constant appropriate 
to six dimensional AdS spaces, n 2 = 87tGq, and the induced metric on the 4-brane is g ab = 


S ab is the stress tensor of an infinitely thin brane located at the cutting surface. The stress- 
tensor so defined may be obtained from the Israel matching conditions 
calculation yields, 


25 . A straightforward 



See 


-2 


H 


H 


-2 


+ 7 ) ^ 

( 12 ) 

\a'\ 

4— gee- 
a / 

(13) 


Here and in the following uppercase latin indices indicate six dimensional coordinates, while 
lower case greek indices specify the coordinates parameterizing the directions along the 3-brane. 
Lower-case latin indices similarly label directions parallel to the 4-brane i.e., x a = ( x^,9 ). 


A crucial point for this model is that the extra term f'/f in eq. (| 12 |) relative to (h3|) , though 
small, is nonvanishing, and therefore it is impossible to interpret the 4-brane stress tensor as 
being due to a pure tension. Were we to do this, thus making S 0 ° = S g e , the 4-brane would be 
forced to go to r = 00, and we would lose the compactihcation of the extra dimensions and the 
localization of gravity. 


There are several kinds of physics on the 4-brane which would naturally involve the required 
difference in S'g 0 — S g e . One is to imagine that the gluing surface is composed of multiple branes. 
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As discussed in ref. 0, one could consider the superposition of the stress-energy tensors of a 
four-brane wrapping the internal circle and a three-brane which is smeared over the internal 
circle. Indeed eqs. (^,|1^) then take the form [0, 



See 



Ta See 




(14) 


where Lg is the proper period of the circle parameterized by 9 at r = R. 

Another very physical possibility is that the difference between Tq and Tg is due to the 
Casimir energy of any massless fields which are confined to the 4-brane (H|. For these the 
stress-energy tensor will take the form 


- hi + f|) v- 


s,„ = 


Ti ■ I) 9W 


(15) 


with some dimensionless coefficients Co and eg. To the extent that the trace anomaly vanishes 
(which is the case at one loop, since the 4-brane is odd-dimensional), the Casimir energy satisfies 
the condition g ab S a b = 0, which implies Cg = 4 c 0 . 

In a static background, either of these stress-energy tensors are trivially conserved on the 
4-brane. But when we discuss dynamical perturbations of the static space, conservation of 
stress-energy will yield a nontrivial constraint on the components of S M „. This will be discussed 
in section 3.3, where we show that the ground state of the radion modes is tachyonic for the 
smeared 3-brane model, but massless for the Casimir model. 

The issue of stabilization is closely related to another potential problem with the above 
models. This concerns the order of magnitude of the difference S 0 ° — S g e , which is required in 
order to acheive a(R) ~ 10 21 as is needed to solve the hierarchy problem. This requires 


c o 


c 0 


Sn° + SJ 


a(R ) 


-5/2 


10 


-53 


(16) 


which appears to be extremely fine-tuned. From this standpoint, only the Casimir effect can 
be considered to be natural, since its dependence scales precisely like a(R )~ 5//2 . However 
the Goldberger-Wise stabilizing held makes it unnecessary to have nonvanishing S 0 ° — S e e , as 
was shown by |L4|: with the scalar it becomes possible to achieve an exponential hierarchy 
even when S 0 ° = S g e . It is interesting that we are able to both determine the size of the 
extra dimension and stabilize the radion using the same scalar held. In the 5D RSI model, the 
two phenomena are necessarily tied together, but not so in 6D. The fact that the size of the 


extra dimension is determined by the value of S 0 


c e 


does not prevent the instability we 


will demonstrate, so it is not obvious that introducing a new effect to determine the size of R 
should stabilize the system. Nevertheless we shall show that it is true. 
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3 Stability Analysis 


In the original model of Randall and Sundrum with two branes, fluctuations of the metric were 
decomposable into Kaluza-Klcin modes. Most notably the spectrum included a zero mode 
which was bound to the brane and served as the graviton in a four dimensional world. The 
remaining excitations formed a tower of massive modes which were fully five dimensional and 
had very little support near the brane. In the AdS soliton model presented here we will find 
a very similar story emerging with a few differences. As in RSI, the spacetime constructed in 
the previous section is finite and one can view the graviton fluctuations as linearized gravity in 
a box. This implies that the spectrum of gravitons will again be discrete. Another difference 
from RSI is that the fluctuations of the metric are now more complicated owing to the greater 
complexity of the background metric. With only a single extra dimension, the only degree of 
freedom for the radion mode is the distance between the two branes, since any apparent ripples 
in the dr 2 metric component can be gauged away by a coordinate transformation. With two 
or more extra dimensions this is no longer the case, and the radion too has a KK tower of 
excitations. 

By virtue of the symmetries of the geometry at least four of the metric modes must be exacly 
massless. First, there are two massless states which correspond to the massless 4-D graviton 
which is ensured by the model’s unbroken Lorentz invariance in the directions parallel to the 
3-brane. Second there are two states making up the massless 4-D spin-one particle, which is a 
consequence of the isometry 6 —> 6 + c of the extra dimensions. 

The counting of massless modes may be further sharpened as follows. If gravity is indeed 
localized on the 4-brane, we would expect to find a total of five zero modes appropriate to the 
five independent fluctuations of a massless spin-two particle in 4 + 1 dimensions. Since we will 
find below that the radion generically has a nonvanishing mass in this theory (either tachyonic 
or real by adding the appropriate scalar), there are in fact only four zero modes bound to the 
brane. These will have a natural interpretation, at energies below the mass of the radion, as a 
3 + 1 dimensional graviton and a 3 + 1 dimensional massless vector field. 

The analysis will proceed by linearizing eq. ([□]) around the background of the AdS soliton. 
In particular we will consider fluctuations of the six dimensional metric which are given by 
Qmn —* 9mn + h\iN- A feature here is the fact that h MN is a tensor and there is thus a variety 
of polarizations, or graviton modes, that need to be considered. Following refs. |2 
can divide the various polarizations of the six dimensional graviton into three categories, (i) 
Transverse traceless polarizations. These are modes which are polarized in directions parallel 
to the Lorentz invariant hypersurface spanned by the coordinates x l . (ii) Vector polarizations. 
These are gravitons whose polarization is of the form e,/), i.e., modes polarized along the circle 
and in the flat piece of the brane. (iii) Scalar Polarizations. These are modes which are 
diagonal but not traceless. It is this last mode which is related to the radion field. 

For the cases (i) and (ii) above we may write the metric fluctuations as Hmn = H&iN(r)e lk ' x 
where H MN {r) is the radial profile tensor and is a 4-dimensional momentum vector with 
k 2 = rj^k^k^ = —M 2 . Further there are ambiguities in the metric perturbations arising 
from diffeomorphism invariance, which we (partially) handle by imposing a “transverse gauge:” 
Hm^ = 0. For massive excitations we may always, via the appropriate Lorentz boost, choose 
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to work in the rest frarnef] so that the momentum can be written as In this case, the 

transversality condition becomes, 

H aiJ k' J = 0 =* H at = 0 Va (17) 

Our implicit notation for the 3 + 1-dimensional Minkowski space coordinates is x ** = ( t,x l ) 
with i — 1,..., 3. 

We do not consider here the Kaluza-Klein modet|] corresponding to angular excitations, i.e., 
around the large extra dimension. One might at first have thought that these had a mass gap 
of order 10 eV since they are modes which are localized on the 4-brane (assuming they are 
radially unexcited) and the circumference of the 4-brane is of order Lg ~ (10 eV) -1 . However, 
if we imagine integrating out only the radial dimension to obtain the effective theory of these 
modes, we find that the fundamental scale is no longer TeV, but rather \Ja(R) TeV ~ 10 13 
GeV, because of the effect of warping. (The kinetic term of these excitations gets the same 
exponential factors as does the angular gradient term: C ~ (0) 2 3 /a + (dg(p) 2 /b.) This effect 
of warping was alluded to in section 2.2. From this point of view, the size of the compact 
dimension looks like (TeV) -1 , whose smallness compared to ~ 10 13 GeV is how the largeness 
of the extra dimension is manifested. We will leave aside these angular KK modes and instead 
consider the radial excitations. To determine the spectrum, we must solve eq. (ED with the 
ansatz (0 as an eigenvalue problem for the mass M. 

The metric fluctuations for case (iii) are considerably more complicated and will be dealt 
with separately in section 3.3. 


3.1 Transverse Traceless Modes 

As explained above, these modes are polarized parallel to the Lorentz invariant directions on 
the brane and correpond to 


HgM — HrM — 0 — Hm V M , (18) 

where the last equality is a restatement of eq. & 

A consistent solution to eq. (ED linearized around the AdS soliton is provided by the 
following ansatz, 

H M n = £ MN a(r)H(r) (19) 

where e MN satisfies the conditions in eq. ([H|) and a(r) is the metric function appearing in the 
static solution above. 

Solving the equations of motion, which come from linearizing eqn. © around our ansatz, 
imposes that the polarization, must also be traceless, 

rF'e tu , = 0 . 


2 Of course when searching for zero modes we are constrained to work with a null momentum vector. 

3 For a discussion of these KK modes as they relate to the stability of the AdS soliton see ref. (23). 


( 20 ) 




Thus eq. (|T8|) describes five independent modes, which can be described as three off-diagonal 
polarizations, e.g., 

£12 = £21 = 1 , otherwise e MN = 0 (21) 

and two traceless diagonal polarizations, e.g., 

Eh = -£22 = 1 , otherwise e MN = 0 . (22) 


For all of these independent polarizations, the radial profile H{r) satisfies the same differential 
equation. Substituting the above ansatze, (|T8|) and (|H|), into eq. ([TT]) and linearizing around 
the AdS soliton background one finds that the radial profile must satisfy 


<PH(r) ( ,a'(r) fWt = n 

dr 2 ^ a(r ) /(r) J dr /(r)a(r ) 2 


(23) 


where primes denote differentiation with respect to r. Here the h-function coming from the 
right hand side has been canceled exactly by similar terms on the left hand side. One should 
note that this is exactly the equation for the transverse traceless modes originally obtained in 
ref. 1[23|. |26H — see also ref. (§[]. Again, this is independent of the form we choose for the stress 


tensor, since by exciting these transverse traceless modes we are not perturbing the size of the 
circle. It is interesting to note that this is precisely the equation describing the propagation of 
a minimally coupled massless scalar on the AdS soliton background 


Here we will determine the eigenvalues numerically using a shooting technique (see ref. ||27|| ). 
For the purposes of numerical calculation we will henceforth restrict to r < R and replace the 
brane by effective boundary conditions on the gravitons at r = R. Obtaining the correct 
solution to this problem is equivalent to determining the correct boundary conditions that the 
radial profile H(r ) must obey at r = p and r = R. At the brane the absence of (^-functions in 
eq. (p3|) implies that H(r) and its first derivative are continuous so requiring an even function 
of r gives the boundary condition H'{r ) = 0. The boundary at r = p is more subtle since the 
metric function /(r) vanishes there, i.e., f(r = p) = 0. This is reflected in the fact that this 
point is a regular singular point of eq. (|23|). So requiring that H{r ) be regular at this boundary 
gives the condition 


dH{r) 

dr 


L 2 M 2 
5 


h MU 


r=p 


(24) 


r=p 


It is now straightforward to see that this polarization of the six dimensional graviton indeed 
has a zero mode. This can be done by setting M 2 = 0 in eq. (|23|) and integrating directly. 
After the first integration we find 


I — const, x a(r) 3 /(r) 1 (25) 

dr 

which can be seen by inspection to violate the boundary condition at the brane unless the 
integration constant is forced to vanish. Performing the second integration just leaves the 
constant solution obeying the above boundary conditions. So we see, referring back to the 
ansatz in eq. (B. that the physical zero mode for these polarizations is 


h-MN — £MN^(r ) 


(26) 
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where e MN obeys the conditions in eq. flT9| ) and we have used our freedom to perforin one overall 
rescaling of the solution to set H{r) = 1. 

In order to obtain the spectrum of nonzero modes we use the numerical shooting technique 
with the above boundary conditions. The mass eigenvalues are a function of the relative size 
of the extra dimension R/p, but in the limit that R/p becomes exponentially large, as desired 
to solve the hierarchy problem, the masses quickly approach their asymptotic values. We give 
these values for the first few KK modes in the following table. We emphasize that there are no 
modes with M 2 < 0 and hence no instabilities in this sector. 


Mode Number 

M' 2 L 2 

0 

0 

1 

16.494 

2 

44.731 

3 

85.545 

4 

138.92 

5 

204.85 


Table 1. Mass squared of the radial graviton KK modes, in units of the AdS curvature radius, in the 
limit of large warp factor. 


Unlike the zero mode, which is localized on the 4-brane, the KK modes are peaked at 
the TeV brane, a phenomenon which is familiar from the 5D RS model. This behavior 
is illustrated for the first three modes in figure 1, where the wave functions are plotted. 



Figure 1. Wave functions for the first three radial KK modes of the graviton. 

3.2 Vector Modes 

The next set of polarizations comes from the same ansatz as in eq. (|l^); however in this case 
the polarization tensor is such that it has one leg in the lorentz invariant directions of the 
four-brane and another on the circle. It has the form 

Eq^ = e = Vfi with k ■ v = 0 and v ■ v = 1 (27) 
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Polarizations of this form contain three independent modes. Substitution into the equation of 
motion (|TT]) and linearizing as before yields 


d 2 H(r ) a'(r ) dH{r) 


dr 2 


+ 


M 2 


a(r) dr /(r)a(r) 


;H(r) = - R)H(r) 


(28) 


For these modes there is a net contribution from the 5-function source term on the right hand 
side of eq. (|ll|). This can be understood from the fact that the metric perturbation we are 
considering is a fluctuation in the gg x components of the metric. The corresponding variation 
of the stress-energy tensor only involves the pieces proportional to T 4 and not those proportional 
to T 3 . In other words for these modes we have effectively that 5T 4 = 5T 0 under the fluctuations 
condsidered in this section. In our formalism this will manifest itself as a nontrivial boundary 
condition at r — R, 


dH(r) 


dr 


r=R 


f ( r ) 

f( r ) 


H(r) 


r=R 


(29) 


while enforcing regularity at the singular point r = p requires that H(p) — 0. As with the 
transverse traceless modes of the previous section, we can obtain an analytic solution for the 
zero mode of this equation by setting M 2 = 0 and performing the integration directly. We find 
an exact solution of the form 

H(r) = ~~ + b (30) 

where c and b are arbitrary (dimensionful) constants of integration. Imposing the boundary 
condition in eq. (|29|) gives 



(31) 


and from this it is straightforward to see that the zero mode is given by 


//(r) = v /(r) (32) 

Fortuitously this solution also satisfies the boundary condition at r = p for all values of c since 
this is precisely where f(r) = 0. Choosing to normalize the wavefunction so that H(R ) = 1 
amounts to choosing c = 5jjf(Rff 1 . Returning to the ansatz in eq. (|27|) we see that the 
physical zero mode takes the simple form 


H fj, T 


f ( r ) 

f(R) 


(33) 


which is indeed peaked at r — R. In order to analyze the nonzero modes we again turn to 
numerics and find a positive definite spectrum which implies that no instabilities are caused 
by exciting these vector modes. In the table below we again present, in the limit of large warp 
factor, the first few eigenvalues in this spectrum. 

Like the spin-2 modes, the vector KK modes are also peaked near the TeV brane, although 
their wave function vanishes precisely there. The first three modes are shown in figure 2. 
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Mode Number 

M 2 L 2 

0 

0 

1 

25.001 

2 

59.752 

3 

106.91 

4 

166.61 

5 

238.80 


Table 2. Mass squared of the radial graviphoton KK modes, in units of the AdS curvature radius, in 
the limit of large warp factor. 



Figure 2. Wave functions for the first three radial KK modes of the graviphoton. 

3.3 Scalar Modes 

Here we will consider modes which fall into category (iii) above. There is in fact only a single 
mode (at the lowest KK level), although it corresponds to simultaneous fluctuations of different 
components of the metric. Most importantly this mode includes fluctuations of both goo and 
g rr and thus couples fluctuations of the radial size to fluctuations in the size of the compact 
circle. For the radion, it is convenient to switch to the polar coordinates introduced in eq. (0), 
where the equations take a somewhat simpler form. We start by writing the perturbed metric 
as 

ds 2 = a(r, t)g^ u dx^dx v + b(r, t)d0 2 + c(f , t)df 2 (34) 

where the perturbation around the static background corresponds to 

a(f,t) = e -M r )-M r t) = ao{r)e- Al W 

b(r,t) = e-BoW-Bi-m = bo^e- 31 ^ 

c(r,t ) = l + Ci(f,t) (35) 

and we have now included a subscript ‘ 0 ’ to denote the metric functions of the background 
around which we are expanding, i.e., the AdS soliton. Here we have gone to the 4-D rest frame 
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of the fluctuation 


A 1 (r,t)=A 1 (f) Re(e~ iavt ) (36) 

and similarly for F\ and G\. This is generically valid since for arbitrary values of the parameters 
there will be no zero mode. We will find however that there is a special case in which there is a 
zero mode and hence no rest frame. For this case we have checked that the procedure outlined 
here gives the correct result. 


We can unify the models of the 4-brane stress energy at f 
section 2.3 by writing 


S, 


flV 


s. 


ee 



9 fj.u 


W g^v 


gee = Ve gee 


R which were discussed in 


(37) 


where Lg — f dO\/b is the circumference of the compact 4-brane dimension. The smeared 
3-brane model has a = 1 and Tg = 0, while the Casimir effect model has a = 5 and T 3 , Eg A 0. 
Notice that the metric tensor elements appear within the definition of the circumference Lg. 
Also T 3 really has the dimensions of a 3-brane tension only when a = 1. 


We can now write the Einstein equations for this metric ansatz. Since it will later be 
necessary to add a bulk minimally-coupled scalar held, we include it here, although at first 
we shall carry out the analysis with no scalar. Ignoring terms like a 2 which would be higher 
order in the perturbations, the (ft), (ft) + {ii), (rr), {98) and (tr) components of the Einstein 
equations are 


3 a" 3 a' b 1 3 a' c 1 lb" 1 / b' \ 2 1 b'd 

2 a 4 a b 4o c ^ 2 6 4^6/ 46c”*~ 


- -k 2 (c[ A + V{<!>)] + \<$ 2 + E 0 V~c5{r - R) + T b ^-8{f)^j (38) 


2 - + - + - 

a b c 

3 fa'\ 2 + a'b' 1 9b ^d\ 
2 \a J a b 2a 3 \b a) 

2 A + i fEV _ Ed _ J_ f 3 4 + 

a 2 \ a J a c 2ao \ a cj 

a' a a' b a 1 c a' b' bb' c b' 

6 6 ----3— + 2------- 

a a a b a c a b b b c b 


0 

(39) 

-K 2 (c [A + V{(j))] - i0 /2 ) 

(40) 

— k 2 (c[A + V (</>)] + + V$-\/cS{f 

-R)) 


(41) 

4K 2 <j)(j)' 

(42) 


where T b is the tension of the 3-brane at r = 0 and primes denote d?. 


The next step is to linearize the held equations in the dynamical perturbations. We use the 
relation Ai = —m 2 Ai, and similarly for the other perturbations, where m r is the sought-for 
radion mass. Expanding the (ft) + {ii), {tr) and (rr) Einstein equations, respectively, to hrst 
order gives 

m 2 {2Ai T B\ — Ci) = 0 (43) 
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- -A'qCx) + 


B' 1 -A' 1 --(B' 0 -A' 0 )(B 1 + C 1 ) 


— 2k 2 0q0! = 0 


(44) 


(3^o + Bq)^ + A' 0 B[ - A' q + B'^j 


ml 


C 1 --^(3A 1 + B 1 ) 
2 a 0 


K 

T 


BV 

i>oCi + 2—0! - 20o0! 


= 0 


and similarly for the scalar field we find 


0" - (2A' + B')0' - -(4Ai + B[ + C')0(, - 


' dV d 2 V m r 

, Cl dj + W * ~ tW 


(45) 


- -54,, = 0. 


(46) 


We note that even for a massless mode, the energy is nonvanishing, so that eq. (43) would 
still provide a constraint among the components of the perturbation (we would then have to 
consider its spatial momentum too). Moreover (f|4|) is the time-integrated form of the (tr) 
equation, which we have organized in a form whose purpose will become clear momentarily. 
We have not written the (tt) nor (99) components since these are not independent equations: 
(tt) can be obtained from (tr)' and (tt) + (ii), which follows from the Bianchi identities; and 
(99) follows from combining (tr), (rr), (rr)' and 0 equations in a manner which is not obvious, 
but which could be anticipated, since first order constraint equations, like (rr) and (tr), must 
be consistent with the second-order dynamical equations, (tt), (ii) and (99). 

For the remainder of this section, we will assume there is no bulk scalar held present. Effects 
of the bulk scalar will be considered in the next section. 

The boundary conditions at the 4-brane come from integrating over the delta functions 
to find the discontinuity in the first derivatives, and using Z -2 symmetry across the brane to 
identify, e.g., Aa'/a = — 2a'/a at r — R. At zeroth order in the perturbations, this gives the 
jump conditions (TT2|.p~3|) for the static solution. Expanding to first order, and assuming that T 3 
is a constant (more about this below) we find 


A = 


B[ - A[ 



(B'„ - A)Bi 


= ;« - A) (Ci + cB,) 


(47) 

(48) 


In addition, we must consider the boundary condition at the 3-brane at the center, f = 0. With 
two extra dimensions, the effect of such a point-like source is to introduce a conical singularity 
and a corresponding deficit angle, as we have discussed in section 2 . 1 . Since the defect is 
unchanged by perturbations around the static solution, we insist that the deficit angle does not 
vary. If we consider a circle with f — e around the 3-brane, with circumference L and physical 
radius D, we therefore demand that L/D he invariant in the limit that e —> 0: 


lim 5 (— 

6-0 \DJ 


- -4^) 

V JoVE / 


2D 


(fli + Ci)!^ 


f =0 > 


(49) 
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in other words, Bi + C\ = 0 at f = 0. In addition, we expect A\ to vanish at r = 0. The latter 
can be shown to be satisfied from B\ + C\ — 0 combined with the bulk equations of motion, so 
it gives no additional constraint. 


We must pause to discuss an important point, concerning the counting of boundary condi¬ 
tions versus independent differential equations. Using the algebraic constraint ([43]) to eliminate 
Ci, we have two first order o.d.e.’s for A\ and B\. On the face of it, our system looks overcon¬ 
strained: there are three boundary conditions! But in fact the system is not overconstrained. 
Rather, there is a constraint on the stress-energy tensor on the 4-brane, due to its conservation. 
By computing S 0a a = 0, where S a b is the surface stress energy on the 4-brane, we obtain 


dTz 

dt 


B\ 


(U(l 


a) + n) 


Integrating this, we see that unless the constraint 


(50) 


n = (a - m 


(51) 


is satisfied, then T 3 must have had extra hidden dependence on B (hence the circumference 
of the circle Lg ) beyond that which was explicitly assumed. If T 3 is truly constant, then 
any physical model of the stress-energy must satisfy (0). This is true for the model which 
corresponds to delocalizing (smearing) a 3-brane around the circular dimension of the 4-brane, 
since there a = 1 and T 3 = 0. And it tells us that the Casimir energy model with a = 5 must 
satisfy T 3 = 4T 3 , as we saw earlier was indeed the case for massless particles. With any such 
choice, it is easy to see that there are not really two boundary conditions at the 4-brane; rather, 
imposing the first b.c. (|47|) , together with the ( tr ) equation (]44|) , implies the second b.c., (|48|) . 
The result of this discussion is that it suffices to impose just one b.c. at r = R, say (|47D, which 


using (51) can be written more simply as 


A = 2^A 


OL 


1) 


(Bi - A' 0 )B,. 


(52) 


We solved the above system of equations for the radion numerically, for the case of a — 1 
(the general result for arbitrary values of a will be given below) and we find that it has a 
negative value of m 2 r —it is a tachyon. The graph of its dependence on the ratio of warp factors 
between the 4-brane and the 3-brane is shown in figure 3a. If we normalize a(0) = 1 at the 
3-brane, then in the regime where the hierarchy is large (right hand side of the graph), the 
radion mass squared depends on a[R ) (a evaluated at the 4-brane) as 

m 2 r = — 20 L~ 2 a(R)~ 3l/2 (a = 1 case only), (53) 

where we recall that L is the AdS curvature radius. Since L 1/TeV and a{Rf' 2 2* 10 32 to 
solve the hierarchy problem, we obtain 


\J—m? r ~ 10 3 eV. 


(54) 


This is well above the present Hubble scale, so we would have noticed the expansion or con¬ 
traction of the extra dimension due to the change in the strength of gravity, if eq. (|54l) were 
true. 
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Figure 3. (a) Log of minus the radion mass squared versus log of the warp factor, for the case a = 1. 
(b) Outer (solid) lines: the radion wave function for a large value of the warp factor; inner (dashed) 
lines: wave function of the first KK excitation of the radion. 


The radion wave function in the large warp factor regime is shown in figure 3b. If the 
4-brane is moved even farther away (larger warp factor), the wave function retains the same 
form, since it stays flat in the region of large f. Although these plots were made for the case 
a = 1, we find that the wave function looks essentially the same for all values of a. We see 
from its functional form that near the 3-brane Ai = —B i, and since C\ is constrained to be 
3Ai + B\ , therefore C\ = —2 B x in this region. For larger values of r we have A 1 = 1/3 and 
Bi = —2/3, so C i is extremely small throughout most of the bulk. Nevertheless its integral is 
nonvanishing, so the radial size of the extra dimension, which is given by f ( fdrC(r), changes 
in response to the instability, and it does so in the same sense as the size of the compact 
dimension, because of our choice of signs in the defintions of Bi and C\. That is, the instability 
is a simultaneous growth or shrinking of the radius together with the circumference. Either 
direction is a possibility, since the static solution is analogous to sitting on the top of a hill: the 
ball can roll down in either direction. The situation is illustrated in figure 4. Also shown there 
is the fact that the relative size of the brane directions, x grow or shrink in the opposite sense 
relative to the extra dimensions. In the case where the extra dimensions grow, the endpoint 
must be the AdS soliton solution with no 4-brane, since this has been demonstrated to be the 
minimum energy solution which is stable N23|. In the case where they shrink, the 2-D surface 


presumably degenerates into a point. 

In comparing the radion in this model to that of the 5-D Randall-Sundrum model, we can 
notice several similarities and differences. Similarly to the 5-D model P8L |29| , in 6-D the radion 


is an admixture of the radial and brane metric components, such that oscillations of the radial 
size are accompanied by fluctuations in the scale factor of the 4-D universe which are 180° out 
of phase. But in 5-D, the radion was exactly massless in the absence of stabilization, whereas 
in 6-D it is a tunable parameter. Another difference is that, whereas in 5-D the radion has no 
tower of KK excitations, in 6-D it does. The mass gap is of order 1/L, i.e., the TeV scale. The 
first few eigenvalues, for large values of the warp factor, are given in table 3. One notices that 
these masses are systematically smaller than those of the graviton and vector modes. Thus the 
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Figure 4. Illustration of the instability. A given static solution, shown in the center, is unstable 
toward growth (left) or shrinkage (right) of the extra dimensions (funnel), accompanied by the opposite 
behavior of the directions x ^ within the 3-brane (shown as a sphere). 


Mode Number 

m 2 r L 2 

1 

1.0188 

2 

6.1512 

3 

12.748 

4 

21.311 

5 

44.437 

6 

75.588 


Table 3. Mass squared of the radion KK modes, in units of the AdS curvature radius, in the limit of 
large warp factor. 

radion excitations would be the first signs of new physics from this model (once we have made 
it viable by stabilizing the radion) to appear in accelerator experiments. The wave functions 
of the second and third excited states are shown in figure 5, while that of the first excited state 
appears in figure 3b. Similarly to the spin 2 excitations, the modes above the ground state have 
wave functions which are exponentially strongly peaked on the TeV brane. The wave function 
of the radion ground state, on the other hand is relatively flat throughout the bulk. This plays 
an important role in its couplings, as we will discuss in the next section. 

We can understand the preceding results for the radion ground state mass and wave function 
analytically, and generalize them to arbitrary values of the 4-brane stress-energy parameter a. 
Toward this end, we first convert the coupled first order equations into a single second order 
equation. The form of the equations suggests that a natural dependent variable to consider is 
the linear combination H = 2>A\ + B \. The variables A\, B\ and C\ can be expressed in terms 
of H using eqs. ([43|) and (f44|), where for later convenience we continue to show the effect of a 
bulk scalar held, even though we set it to zero for the present: 
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Figure 5. Solid lines: wave function of second excited state of the radion; dashed lines: the third 
excited state. 

1 


Ai = 

B 1 = 
C x = 


2 B' 0 L 


-H' + (A' 0 + B' 0 )H- 2k 2 0o0i 


2 B' 0 
1 

2K L 


H'- A’+-B')H + 2k 


.2 jJ 


2 0 <dl 


H' + (£' - A')# - 2k 2 ^fa 


(55) 


Substituting these into the remaining field equation (|45P , we obtain 


H" + 


QA'n - 3 B' n 2 - 6AqB' 0 - 2k 


,2u2 


2 B’ n 


-H' + 


{Bq — A' q )(?>A I q + 2B' 0 A' 0 — k 2 0q 2 ) , m 


B' 


+ ^)H 

a 0 I 


~ 4 « 2 (w " k 2 k v ) h = 0 (56) 


and the boundary conditions at r = R and r = 0, respectively are 

H’ = 


A' 0 A-^B' 0+ 


2 A' 2 
0 


1 


A',) + -Bn ) ( A'n + 


3 (B’ 0 - Aft 
1 - 


a 


B’ 


2B' 0 ml 


2 2 i / 2 

3 + 3(^-^ 0 )a 0/ 




r = R 
r = 0 


H' = A' 0 H — 0, f = 0 (57) 

The radion mass squared comes into the boundary condition because, in the process of elimi¬ 
nating Ai, A \, and B[ in favor of H' and H\ it is necessary to use the bulk equation (|45|), 
evaluated at the 4-brane. 

Now it happens that the important features of the radion ground state can be deduced from 
approximating the above equations by the form which they take in the asymptotic region of 
r ~ R when R is large. In this region, we have 


A' 0 = B' 0 = — k- B' 0 — A' 0 = -8 ke~ 2kf = 5 A' 

5 


( 58 ) 
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and the equations simplify considerably: 


H" - -A'H' + 15^6A' + — 

2 l a 0 , 


H = 0 in the bulk; 


(59) 


2-'2 3 


mt 


=A + g( 5 ~a)M'|J7'=[£ + 


-AAu A '\H at f = fi, (60) 


where now we take A' to have the constant value —4k/5. The terms 5A' and m 2 /aQ are both of 
order e~ 2kr in the large f region (as we will verify self-consistently), so that in the bulk equation 
(|59|) they can be ignored compared to the other terms. The solution in the bulk has the form 

H(f) = ci + c 2 e- (6fc/5)r ~ + SH, (61) 


where 5H represents the small effect of the parenthetical terms we have ignored. The latter give 
rise to a negligible effect on the bulk solution, 8H -C H. However, the small terms proportional 
to 8A' and m 2 /ao cannot be neglected when applying the b.c. at the 4-brane. In fact, this 
equation, (|60D, can be used to solve for the radion mass, which in the absence of the scalar field 
gives 


m 


2 

r 


5 — a , f 3H' / 

\2H- 2A , 


f=R 


- (° - 5 > J|y" 2<r< “ /5)S 


(a - 5)-W 3,2 (R) + O(a 0 (R)SA r2 ) 

ol- 5 A / A \ 

2 ^ /2 («) WA)J 


(62) 


The final expression assumes that ao is normalized to unity at the TeV brane, and uses the re¬ 
lation (|H) between A and L. (The intermediate factor of 2 4 / 5 comes from ao(R) = cosh 4 ^ 5 (kR).) 
Interestingly, the value for rri 2 which we so obtain is completely insensitive to the details of 
tL e -( 6 fc/ 5 )R, [ ess since all of these are much smaller than A'. We are therefore able 

to give a very accurate analytic estimate for the radion mass squared, when the warp factor is 
large. The small magnitude of rr? T is seen to be a direct consequence of the value of B' 0 — A' 0 
in the static solution. This expression agrees with our previous numerical results for a = 1 
(and we have also checked it numerically for other values of a). In the limit that the 4-brane 
goes to infinity, so that the full AdS soliton is recovered, the radion becomes massless, but is 
not normalizable. Thus it does not contradict the fact that the uncut AdS soliton is a stable 
solution. 


Interestingly, the radion mass vanishes almost exactly in the case where the anistropy of 
the 4-brane stress tensor is provided by the Casimir energy and pressure of fields living on 
the compact extra dimension. In this case the relevant energy density scales like L / 5 , i.e., 
a = 5, as expected from dimensional analysis. This is the unique case where no dimensionful 
parameter is introduced in the anistropic part of the 4-brane stress tensor, which is the part 
that also controls the position of the 4-brane, and hence how large the extra dimensions are. 
Curiously, the mass does not vanish exactly when a = 5 because of the 0(8A' 2 ) correction, 
whose coefficient turns out to be (11 + ct)/ 8 . However, the natural size of this contribution 
to m r is of order 10 ” 10 eV, which is far below experimental limits on scalar-tensor theories of 
gravity. 
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4 Radion stabilization and phenomenology 


In this section we show how to increase the mass of the radion through using a bulk scalar field, 
and discuss the implications of the model for collider experiments, tests of the gravitational 
force, and cosmology. 


4.1 Stabilization by a bulk scalar field 

We have found that the radion can be massive, massless, or tachyonic, depending on the value 
a which controls the dependence of the 4-brane stress-energy on its circumference. In the latter 
two cases (a < 5), it is certainly necessary to increase the radion mass squared so that we have 
a stable universe, with Einstein gravity rather than scalar-tensor gravity at low energies. In the 
5-D RSI model, this was achieved by Goldberger and Wise [0 by adding a bulk scalar field, 
whose VEV’s at the two branes were constrained by potentials on the branes to take certain 
values. The bulk scalar then acts like a spring between the branes, whose gradient energy 
becomes repulsive if the branes get too close, and whose potential energy (from m 2 (f) 2 ) causes 
attraction if the branes separate too much. We expect that the same mechanism should work 
in 6-D. 


Scalar fields in AdS have solutions which are exponentially growing or decaying toward the 


ultraviolet cutoff brane. Ref. 114] studied these solutions and found the approximate behavior 


0 (r) = <fi+e a+r + 0_e° 


(63) 


where a± = —k ± \/k 2 + m 2 . Near the 3-brane, where the space does not look like AdS, 
the behavior is different; 0(f) ~ 0 O (1 + rn 2 f 2 /4), but this will not be very important for 
understanding the effect of the scalar since most of the volume of the extra dimensions is near 
the 4-brane. For generic boundary conditions, the growing solutions dominate, and it is a 
good approximation to neglect the decaying ones. The main point is that the most natural 
configurations are ones where 0(0) < 0(R). 

Before doing any analytic estimates, we solved the entire system of Einstein equations 
numerically, to find the effect of the scalar held on the radion mass. There are three kinds 
of corrections to consider. First, the scalar held induces a small back-reaction on the static 
solutions, A 0 , B 0} determined by the zeroth order truncation of the Einstein equations ( fk^|42l ). 
This effect has been analytically computed in |L4|. Second, the background scalar configuration 
couples to the fluctuations of the metric. This arises solely through the term n 2 (j)' 0 2 Ci of the 
perturbed (rr) component of the Einstein equations, (|45|). We will see that this is the really 


important effect for stabilizing the radion. The third kind of correction is from huctuations of 
the scalar held, which can mix with the radion. These are governed by the perturbed scalar 
held equation (fib|). 


-(4A'+R')0' 1 . 


W- 


1 dV 

25^00 


((B' 0 -A' 0 )H + H r - 2« 2 0'0i) 


< 9 2 t / rr? 

^0 1 + ^0 1 = O (64) 
0 (p z a 0 


where V = \m 2 (jr is the bulk potential. 
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Our numerical results demonstrating the stabilization of the radion are shown in figure 6. 
We considered scalar field configurations with 0 = 0 at the 3-brane and varied the value of 0 
at the 4-brane, showing that rn 2 (in the tachyonic case a = 1) becomes positive for sufficiently 
large values of (f>(R). (Treating 0(i?) as an adjustable parameter can be justified by imagining 
that we have stiff potentials for 0 on the branes, fixing their boundary values to whatever we 
desire.) We checked that these results are quite insensitive to whether the fluctuation of the 
scalar are included. The mixing between the radion and 0i was found to be negligible. 


0.03 


0.02 

2 t 2 

m r L 

o.oi 


0 


Figure 6. Dependence of rn 2 on the value of 0 at the 4-brane. 



The behavior shown here can be easily understood by generalizing our previous derivation 
of the radion mass to include the effect of the scalar. The equation of motion and boundary 
condition for H , in the large-r region, were given in ( |59| - |6C)| ). We can take m 2 of the scalar to 
be small, so that its effect can be neglected in the bulk equation and the approximate solution 
(M3) is still valid. Now when we solve the b.c. for the radion mass, we obtain the previous 
expression plus a new term, 


o 

m r 


5 — a 


6 A' 



HA! 


(65) 


whose origin can be traced to the extra term in the (rr) Einstein equation. The first term 
is also changed by the presence of the scalar held, because of its back-reaction on the static 
metric. However, using the results of ref. [|T|] who computed this back-reaction, we find that 
8A’ is still of order e~ 2kR . Therefore, since is of order e 3fcr / 2 , the new term on the r.h.s. of 
([65]) is the dominant one. The fact that has the correct sign (negative) to insure that the 
radion mass squared is positive is not obvious, but by numerically solving for H(r) we have 
verified that indeed ^fj(R) < 0, as we show in figure 7. We thus find that for large enough 
values of 0(i?), the radion mass is 


m r 


m 2 (j)(R) kR/5 

kM 2 


MeV 


independently of the details of the 4-brane stress-energy. 


( 66 ) 


21 










It is remarkable that the stabilized mass is not of order the TeV scale, as was the case in 5-D 
RSI [T8|, P0fl . The mass squared is suppressed by the fractional power of the warp factor left 
in the product ao^f ~ a 0 ^ 2 . Recalling that the Planck scale hierarchy was set by a^ 2 = 10 32 , 
we see that the stabilized m r is suppressed by the factor 10 16//3 , giving m r ~ 10 MeV. This is 
precisely the same factor by which a TeV mass particle, transported from f = 0 to r — R, falls 
short of the Planck scale, as we noted in section 2.2. Thus the smallness of the radion mass 
seems to be associated with the additional dilution of the strength of gravity that comes from 
the large extra dimension. 



Figure 7. The radion wave function H{r ) for the ground state, showing that H'/H < 0. 


In ref. [T3], one advantage of having a bulk scalar field was already pointed out: because of 
the back-reaction of the scalar on the metric, the jump conditions at the 4-brane can generically 
be satisfied for large values of f?, as desired for achieving a large hierarchy, without much 
sensitivity to the model of stress energy on the 4-brane. In particular, choices like a = 0 (pure 
tension) or a = 1 (smeared 3-brane), which by themselves could not have yielded a satisfactory 
value of R , become viable in the presence of the bulk scalar. The only requirement for getting 
sufficiently large R is that the scalar mass should be somewhat light since 


R 


k 


m 


2 ' 


(67) 


Hence one needs k/m ~ 8 to get the desired hierarchy. 
D realization of Goldberger and Wise 


A similar relation occurs in the 5- 
This requirement is quite compatible with the 
parameters we need for generating the radion mass, as in retrospect one would have expected. 
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4.2 Couplings of the radion and its excitations 

In the 5-D RSI model, the stabilized radion has a TeV scale mass and TeV suppressed couplings 
to standard model matter |H|. Were the couplings of our MeV-scale radion so large, it would 
easily be observable in low-energy experiments and possibly affect the cooling of supernovae. 
Here we show that the couplings are actually Planck scale suppressed. 
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Computing the 4-D effective Lagrangian for the gravitational fluctuations at quadratic order, 
we obtain 

1 


£n — 


3 k 2 


d6 dr ao(r)Jbo(r) H 2 (r, t) + 2R{(r, t) + A 1 (0,t)T fl 


( 68 ) 


where T IJ /t is the trace of the 3-brane stress-energy tensor, representing the standard model, 
H(r,t ) = H(r)<po(t), Bi(r,t) = Bi(r)(p 0 (t), etc. Here tpo(t) represents the 4-D ground state 
radion field, and H(r), Bi(r), Ai(r) are the corresponding wave functions found in the previous 
section. Since they are nearly constant throughout the bulk, we can take them out of the integral 
and perform it to obtain 

C 0 ~ M 2 (fl + voT* (69) 


This shows that the canonically normalized radion field ground state has Planck-suppressed 
couplings to TeV-brane matter. This differs from the behavior of the radion in the 5-D RSI 
model. There, the wave function of the radion is exponentially peaked at the Planck brane, 
which overcomes the exponential warp factor in the measure to give £o rv./ (TeV)VS + <PoT^ 
instead. The flatness of the radion wave function in the present case accounts for its weak 
couplings to the TeV brane. 


The KK excitations of the radion are exponentially peaked on the TeV brane, on the other 
hand. We can understand this from the asymptotic form of the bulk equation of motion (|59|); 
since the mass is no longer negligible, the solutions behave like H{r ) = C 2 e _6fcr / 5 , with the 
constant piece c\ equal to zero. The integrand of (^) behaves like e - 12kr / 5 + 12kr / 5 — 0(1), so 
we obtain 

M 4 rv 

C n ~-±Rvl + V>nT% + (70) 

Hence the coupling of the radion excited state is suppressed only by the small factor ( MqR) 1 ^ 2 ~ 
v^60 relative to the TeV scale. The radial KK gravitons have similar couplings, but larger masses 
(compare Tables 1 and 3), so the radion excitations would be the first signal of new physics in 
collider experiments. Heavy radions could be copiously produced in the s-channcl at the LHC, 
through gluon-gluon fusion events due to the QCD trace anomaly contribution to T 1 '^. 


4.3 Gravity and cosmology 


Although perhaps less physically motivated, models of the 4-brane stress-energy with a > 5 
predict that the radion mass squared will be positive even without a bulk scalar, and that 
its magnitude is in the ICC 3 eV regime. This is within the reach of Cavendish-type tests 

H 


of submillimeter gravity 


and will be even more accessible to upcoming versions of the 


experiment which will have improved sensitivity. 

ft may also be possible to achieve this situation without appealing to exotic forms of matter 
on the 4-brane. The radion mass will get radiative corrections from its couplings to matter. 
Since the radion couples to the trace of the stress energy tensor on either brane, the heaviest 
particles will contribute the most strongly. Considering matter which is on the TeV brane, we 
can estimate the size of the one-loop correction as 


m. 


r,l—loop 


TeV 4 

M 2 

p 


(71) 
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The numerator comes from the fact that the TeV scale is the cutoff on the 3-brane, and the 
heaviest particles will have masses of this order, whereas the denominator is due to the fact 
that the lowest mode of the radion has Planck-suppressed couplings. This argument could be 
upset if the 4-brane has massive particles which are much heavier than the TeV scale, since by 
the same argument these could apparently make the radion very heavy and presumably would 
destabilize the hierarchy which we have achieved. It may be necessary to assume that there are 
only massless particles on the 4-brane to avoid this. 

On the other hand, if we allow heavy particles to exist on the 4-brane, their natural mass 
scale is \Ja(R) TeV rv./ 10 13 GeV. It is interesting that this is the right order of magnitude for 
generating the obsserved primordial density fluctuations from the simplest model of chaotic 
inflation. This is an advantage of the present model over the 5-D RS model, where a ~ 10 13 
GeV would look unnaturally light were it living on the Planck brane, and of course too heavy 
to exist on the TeV brane. 


5 Summary 


We have focused on the simplest and most direct generalization to six dimensions of the 5-D 
Randall-Sundrum two-brane model: the AdS soliton model, with the TeV 3-brane at the center 
of the azimuthally symmetric extra dimensions, and a 4-brane cutting the space off at some 
finite radius. The model has many features in common with its 5-D predecessor: the geometry 
is highly warped and very close to AdS in the region far from the 3-brane, the graviton zero¬ 
mode is localized on the hidden brane, while the radial KK excitations are localized on the TeV 
brane and have a TeV mass gap. In both models, the radion can easily be stabilized by the 
Goldberger-Wise mechanism, using a bulk scalar field. 


However, there are also some quite distinctive differences. The hierarchy between the Planck 
and weak scales, while generated mostly (2/3) by warping, is also partly (1/3) due to the 
exponentially large size of the compact extra dimension |I3|, giving it some features in common 
with the large extra dimension proposal. The mass scale at the 4-brane is not the Planck scale, 
but it is suppressed by the size of the large compact dimension to the 10 13 GeV scale. There 
is a tower of relatively light (~ TeV) KK gravitons corresponding to this large dimension. In 
the absence of stabilization by a scalar field, the 6-D model requires some mildly exotic form of 
stress energy on the hidden brane in order to have a finite volume. The 4-brane stress tensor 
generically depends on the size of the extra dimension as Lg a with some model-dependent 
number a. 


Most of these features were already known; in the present work we computed the spectrum 
of metric perturbations, including the graviton and graviphoton modes, and we found the 
unexpected new result that the radion is not necessarily massless, but has a mass squared 
which depends linearly on a and the negative bulk cosmological constant: m1 ~ (5 — ct)A 
TeV/Mp. Only for the special case of Casimir energy on the 4-brane (a = 5) is it massless. 
For smaller values of a it is tachyonic, and the spacetime is unstable. Its couplings to the TeV 
brane are Planck suppressed rather than TeV suppressed, due to the different behavior of its 
wave function relative to the 5-D case. Once stabilized by a bulk scalar field, the radion mass 
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is not TeV scale, as in 5-D, but rather at the MeV scale. This suppression is related to the 
presence of the large extra dimension which does not feature in 5-D. 

The 6-D model has similar phenomenology to the 5-D model, since the Kaluza-Klcin exci¬ 
tations of the radion behave much like the ground state of the stabilized 5-D radion. However, 
there is a new possibility that the radion is stabilized not by the Goldberger-Wise mechanism, 
but by some form of stress energy on the 4-brane which has a > 5, or perhaps by radiative 
corrections from standard model particles on the TeV brane. In this case the radion mass is 
in the milli-eV range, which is just right for being accessible in experiments which test gravity 
below 1 millimeter. 

The latter possibility would seem to require the absence of massive particles on the 4-brane, 
since radiative effects there should induce much larger corrections to the radion mass. In fact it 
might be necessary to forbid heavy particles on the 4-brane just to maintain the large hierarchy 
we set out to achieve. This is a question which deserves further study. But if it is consistent to 
have heavy fields on the 4-brane, then the fact that their mass is naturally of order 10 13 GeV 
is intriguing for inflation, since this is the right scale for getting density perturbations of order 
10 -5 in chaotic inflation. 

We have left for future work a study of 3-brane fluctuation modes, where the position of the 
3-brane could oscillate with respect to the center of the extra dimensions. These modes, if they 
exist and are sufficiently light or strongly coupled, could be important for the phenomenology 
of this model, since they might induce a coupling between the graviphoton and standard model 
particles. 
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